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On  the  Exponential  Convergence  of  the  h-p  Version  for 
Boundary  Element  Galerkin  Methods  on  Polygons 

1.  Babuska,  B.  Q.  Guo,  E.  P.  Stephan 


1.  Introduction 

Most  research  for  the  boundary  element  method  (BEM)  has  been  carried  out  in  the  framework 
of  the  h-version  [11],  [13],  [28]  where  accuracy  is  achieved  by  decreasing  the  mesh  size  h,  while  keeping 
the  degree  p  of  piecewise  polynomials  fixed.  For  this  method  several  detailed  results,  including 
asymptotic  rates  of  convergence  for  both  first-kind  and  second-kind  integral  equations  are  well-known 
(see  [21],  [28]).  The  basic  idea  of  the  above  convergence  proofs  is  the  observation  that  for  strongly 
elliptic  operators  one  obtains  quasioptimal  convergence  in  the  energy  norm  for  any  Galerkin  scheme 
with  conforming  boundary  elements  (see  [27]).  This  result  has  been  used  in  [25]  and  [26]  to  analyze  the 
p-  and  h-p  versions  with  quasiuniform  mesh  for  BEM  for  some  first  kind  integral  equation  on  a 
polygon  T.  Those  versions  for  BEM  have  been  recently  introduced  (see  [1],  [2],  [18],  [22],  [23]).  In  the 
p-version  a  fixed  mesh  with  constant  h  is  used,  and  accuracy  is  achieved  by  increasing  the  degree  p  of 
the  polynomials  used  as  boundary  elements.  The  h-p  version  combines  the  iwo  approaches.  If  one  uses 
a  quasiuniform  mesh  on  the  polygon  T  and  if  the  singularity  of  the  solution  of  the  integral  equation  is 
located  at  vertices,  the  rate  of  convergence  for  the  p- version  of  the  boundary  element  method  (BEM)  is 
twice  that  of  the  h-version  (see  [25],  [26])  for  some  first  kind  integral  equations.  These  results  have 
been  known  for  the  finite  element  method  where  these  two  extension  processes  (p-  and  h-p  version) 
have  been  thoroughly  investigated  in  a  series  of  papers  [3],  [9],  [10],  [8],  [14j.  Furthermore,  for  the 
finite  element  method  on  a  geometric  mesh  it  has  been  shown  in  [17]  that  under  proper  assumptions 
satisfied  usually  in  practice,  namely  that  the  given  data  are  piecewise  analytic,  the  h-p  version  has  an 
exponential  rate  of  convergence  with  respect  to  the  number  of  degrees  of  freedom  while  the  h-  and  p- 
versions  have  only  a  polynomial  rate.  In  this  paper  we  show  the  corresponding  result  for  the  boundary 
element  method,  i.e.,  if  one  uses  a  geometric  mesh  r"  on  the  boundary  T  which  is  graded  towards  the 
vertices,  then  the  convergence  rate  of  the  Galerkin  solution  for  the  underlying  boundary  integral 
equation  is  exponential.  We  consider  both  the  Dirichlet  a  *  tnc  Neumann  problem  for  the  Laplacian 
in  a  plane  polygonal  domain  and  solve  them  via  first-kinu  ;gral  equations  on  T  with  a  weakly 
singular  and  hypersingular  integral  operator  V  and  D,  respectively.  It  is  known  [13]  that  the  operator 
V  of  the  single  layer  potential  is  a  continuous,  bijective  mapping  from  the  Sobolev  space  H  (T)  onto 
H1/2(T)  if  for  the  conformal  radius  of  the  polygon  T  there  holds  Cap(T)  ^  1.  From  [12]  we  have  that 


the  operator  D  of  the  normal  derivative  of  the  double  layer  potential  maps  H 1 7 2 ( T )  continuously  and 
“1/2 

bijectively  onto  H  (T).  As  one  of  our  main  results  we  show  that  if  the  given  data  are  piecewise 

analytic  then  the  solutions  of  the  integral  equations  are  also  piecewise  analytic,  i.e..  for  given  data  in 

0.1  1,2 

some  countable  normed  space  g  £  Bj  (F)  we  have  u  6  83  (f)  satisfies 


Du  =  (l-K')g  on  r. 


whereas  for  g  6  (T)  D  C°(r)  we  have  €  B^T)  satisfies 

=  (l  +  K)gon  r,  (1.2) 

where  K,K/  are  the  operators  of  the  double  layer  potential  and  its  adjoint. 

k,i 

In  Section  2  we  introduce  the  weighted  Sobolev  spaces  (fl)  (k  >  i  >  0)  ana  the  countable 

l 

normed  spaces  B/?(Q).  We  quote  some  related  trace  and  extension  theorems  from  [6]  together  with  the 

o 

existence  results  of  the  solution  u  €  B/3(Q)  of  the  Neumann  and  Dirichlet  boundary  value  problem, 
respectively. 

In  Section  3  we  consider  the  integral  equation  (1.1)  on  T  governing  the  Neumann  problem  and 

1,2 

prove  the  above  mentioned  existence  and  regularity  result  for  its  solution  u  €  Bj  (T).  Then  we 
present  our  geometric  mesh  on  the  boundary  curve  F  for  an  L-shaped  domain  and  the  corresponding 
conforming  boundary  element  space  Sp(r”)  consisting  of  continuous  piecewise  polynomials  of 
increasing  degree  where  linears  are  used  on  the  smallest  subinterval,  quadratics  on  the  next  larger 
subinterval  and  so  forth.  Then  we  show  (Theorem  3.4)  that  the  Galerkin  solution  up  £  Sp(r")  for 
(1.1)  converges  in  H  (T)  with  exponential  rate  towards  the  exact  solution. 

In  Section  4  we  present  the  corresponding  results  for  the  Dirichlet  problem  by  showing  that  the 
solution  gjj  of  the  integral  equation  (1.2)  belongs  to  Bjj  (T)  and  that  its  Galerkin  solution  \ bp  € 
Sp  '(fa)  converges  exponentially  in  H-1/,J(r). 


2.  Preliminary 


Let  fi  €  R2  be  a  bounded  domain  whose  curvilinear  boundary  dQ  is  a  piecewise  analytic  curve 

M 

T  =  U  I\,  where  T,  is  an  open  arc  connecting  the  vertices  A<  and  A,+1  (Aj  =  Aj^+j).  We  denote  >1 
the  internal  angle  at  A,  by  w,  and  assume  0  <  w,-  <  2x,  1  <  i  <  M. 


Let  Hk(D),  k  >  0  integer,  denote  the  usual  Sobolev  space  (see  [20]),  i.e., 


Availability  Codes 
Avail  and/or 
blst  !  Special 


□  □ 


H‘(n>  =  Wo<S<|i  IID%HLi(n)  =  ll..ll2H,(n)  <  oc) 
where  a  =  (aj.a,).  a,  >  0  integer,  i  =  1,2,  Jet j  =  aj  +  a2,  and 


Dcu 


d'  u  =  u 

dx?'-dx?-  x?ix?2- 


Since  f2  is  a  Lipschitz  domain  any  u  -  Hl(ff)  can  be  extended  to  H*(R2)  (see  (24]).  Therefore 


H*  1A(f)  is  defined  as  the  restriction  on  f  of  functions  in  K*(Q)  for  integral  k  >  1,  i.e.. 


H  '‘(f)  =  {u|r:u  6H*(n)} 


with 


and  for  k  <  0 


l!gl|Ht-1/2(r) 


inf  ||u|| 
u|r=g 


Ek{Q) 


h‘'1/2(D  =  (H  (i"‘/2)(r))'  (dual  space). 

Let  r,(x)  -  dist(x,,AJ.  3  =  (31,3->....3 w)  be  an  M-tuple  of  real  numbers  0  <  3,  <  1.  For 

M  3  -u  k 

any  integer  k  >  0  we  shall  write  3  +  k  =  (^1  +  k,/?2+k,...,^!^+k),  and  $5+l(x)  =  f]  r,  1  (x).  We 

define  the  weighted  Sobolev  space  for  integers  k  and  £,  k  >  £  >  0,  by 

(fi)  =  {u|u  €  H^_1Q  if  £  >  0,  H^>|3  +  |q|_^°uIIl  (Q)  <  oo  for  0  <  £  <  |a|  <  k} 

and  the  countable  normed  space  for  £  >  0 

b{(0)  =  {U  €  H^(fi),  V  k  >  £,  ||-t>5+i_£  D°u||U(fi)  <  Cd‘-£(k-?)! 

for  |a|  =  k  =  £ ,  £  +  1 . with  C  >  1,  d  >  1  independent  of  k) 

The  space  hJ  1/j'^  1/2(r),  (resp.  k,  £  integer,  k  >  £  >  0,  is  the  trace  space  of 

H g  (Q),  (resp.  B$(Q)),  i.e.,  for  any  g  €  Hfl-1/2,^-1/2(r)  (resp.  B^-1/2(r))  there  exists  G  6  Hj  (Q) 


4 


(resp.  B^ffll)  such  that  Gjp  =  g,  and 


k- 1/2, £-1/2 

Hfl  (f) 


inf  1 1 G !  i 
G|r=g 


k.l 

Ha  (fl) 


The  analogous  definition  of  the  weighted  Sobolev  spaces  and  countable  normed  spaces  on  the 
interval  1  =  [a,b],  a.b  G  R.  we  quote  from  [C]: 


For  k  >  £  >  0  integer, 

fc.£  0  .  .  t  *r\  1 

H3  (I)  =  {u|u  €  He_1(I)  if  £  >  0,  ||$-+m_£  U  '(x)||U(I)  <  00  for  0  <  l  <  m  <  k}. 


and  for  l  >  0 


B*(I).-  {u  €  H-  £(I),  V  k  >  £,  !l‘i>i  +  t_^(,)(x)||L_j(I)  <  Cd*"£(k-£)! 

for  k  =  £,£  +  1,...  with  C  >  1.  d  >  1  independent  of  k} 


where  <f>-+t_£(x)  =  fl  *(x).  Mx)  =  lx~al>  Mx)  =  lx~bl.  0  =  (^1^2)-  0  <  3 1-  ^2  <  1- 

1=1  k,£;  ij 

For  any  T,  €  T  the  spaces  Hg  (T,)  and  Bg  (T,)  are  defined  with  the  help  of  a  smooth  map 

k.l  l  k.l  M  ki,  £ 

I  —  T;  via  the  spaces  (I)  and  B^(I).  We  define  Eg  (T)  :=  f]  (T;)  and  B-,(T)  :  = 

M  ^  .....  j=1  . 

n  B y ^. ( T j)  with  l  =  (£1 . £m)  and  0  =  (/?,,/?2,".,i?M),  0 ,  =  (^,,1,^.2).  We  shall  write 

j=l 


3j  >  3 j,  if  3jk  >  0jk,  k  =  1,2,  and  3  >  3  if  3 }  >  3,,  1  <  j  <  M.  For  any  real  number  s  we 
shall  write  0  ■  >  s  (resp.  0  >  s)  if  0 }  k  >  s  (resp.  0j  >  s). 


_ 1/2  £  — 1/2 

It  is  difficult  to  verify  whether  a  function  on  T  belongs  to  the  spaces  Bg  ‘(T)  and 

,  k,l  £ 

Eg  1/2(T).  On  the  contrary  the  spaces  (T)  and  B^(T)  characterize  the  traces  of  functions  of 

k.l  ( 

Bg  (fl)  and  Bfl(fl)  on  T  in  a  precisely  verifiable  manner  (see  Theorem  2.1).  Theorem  2.2  deals  with 
the  extension  of  functions  from  T  to  fl. 


£  *.£ 

Theorem  2.1.  (cf.  Theorems  4.1  and  4.4  of  [6]).  Let  u  G  Bg(fl)  (resp.  (fl)),  £  =  1,2,  k  >  £+1, 

£-1  k,l-l 

then  for  1  <  i  <  M,  u|r  G  B^.  (T,)  (resp.  B'g  (I\))  with 
1  «  *  * 


0 


3,.:  €  (3,  +  ]^-l  k).  )  =  1,2  ./  4  <  J>  +  1  <  1  (2.1a. 

£  t  £ 

oru|p  €  8j  (Tt)  (re.sp.  (T,))  with 

1  1  1  1 

3i.i  e  */ 0  <  J,.J1  +  i  <  l  (2-lb. 

Theorem  2.2.  (cf.  Theorems  4.3  and  4.5  of  [6]).  Let  g  be  defined  o *  T  and  g,  =  g  p  . 

(i)  // g  €  C°(r),  and  g,  £  B£,(I\)  (resp.  k  >  2)  for  0  <  j,  <  s  or  g,  £  B j, t <  T , .  .resp. 

H$  ( r, ) ,  k  >  2)  /or  i  <  3 ,  <  1,  1  <  i  <  M,  then  there  exists  G  £  B  >•.  Q !  •  r?  >p.  H  »•  ( Q i »  such  that 

G|r  =  g- 

(»)  If  g,  e  B^(r,.)  ( resp .  Hjj.  (T,),  k  >  1)  for  0  <  3,  <  or  g,  e  B^'T,  f-sp.  HijT,))  for 
,  1  t.l 

5  <  3,  <  1,  1  <  i  <  M,  then  there  exists  G  £  B3*(Q)  (resp.  H^*(Q))  such  that  (J  j-  =  g. 

In  (:)  and  (it)  3 ’  =  ty»£A  /?*  satisfying 

3'  =  ma x(3,_L2Jn),  3t  J  =  3,  }  -  \  sign  (3,  ;-j),  (2.2) 


for  1  <  i  <  M,  1  <  j  <  2. 

In  the  following  sections  we  assume,  for  sake  of  simplicity,  that  Q  is  a  straight-line  potvgon. 

1.1+ 1 

We  will  make  comments  on  a  curvilinear  polygon  at  the  end  of  the  paper.  By  B3  (T)  we  denote  for 
1  <  i  <  M  the  space  []  B3,(r.)  x  [I 

Q<3i<\  \<3t<\ 

In  this  paper  we  consider  the  Dirichlet  problem  of  the  Laplace  equation 


Au  =  0  in  ft 

nip  =  g 


(2.3) 


and  the  Neumann  problem 


Au  =  0  in  fi 


=  g 


(2.4) 


6 


where  ^  means  the  normal  derivative  with  respect  to  the  unit  outer  normal  n,  and  g  satisfies 
an 


|  gds  =  0 

r 


Combining  Theorem  2.1  of  [4]  and  Theorem  2.2  above  we  have  the  following  theorems. 

Theorem  2.3.  // g  6  B3V)  D  C0(T)  with  8  =  (3xj2,..J^  3,  =  (3,  ..3,  2).  0  <  3,  ;  <  1. 

2 

1  <  i  <  M,  1  <  j  <  2.  then  the  problem  (2.3)  has  a  unique  solution  u  t  B 3 ( Ti )  with  3  given  by 


3,  =  3]  if  3'  >1-5; 


3,  >  1  -  £  </  3’  <  1  -  f- 

"*»  *  l 


where  3  x  satisfies  (2.2) 


Theorem  2.4.  // g  €  B^‘(T)  with  3  =  (3X,32 . &$[)'  K  =  Widish  0  <  h.3  <  J-  1  <  '  <  M- 

2 

1  <  j  <  2.  then  the  problem  (2.4),  (2.5)  has  a  solution  u  €  B5(Q)  which  is  unique  up  to  a  constant 
with  3  given  by  (2.6). 


3.  Boundary  Element  Method  for  the  Neumann  Problem 

We  consider  the  Neumann  problem  (2.4)  with  solvability  condition  (2.5).  As  a  consequence  of 
Theorem  2.1,  we  obtain 

Theorem  3.1.  Let  u  €  B^(fi),  then  €  B^F)  with  3, }  given  by  (2.1). 

A  combination  of  Theorem  2.1  in  [4]  and  Theorem  4.5  in  [6]  leads  to  the  result 

Theorem  3.2.  Let  g  €  B^'V)  satisfying  (2.5)  then  there  exists  u  6  B^(Q)  solving  (2.4)  with  3  9"en 

1,2 

by  (2.6).  Furthermore  u|p  €  B^  (T)  with  3  satisfying  (2.1). 

Proof.  Due  to  Theorem  2.2  there  exists  a  G  €  B^*(D)  such  that  G|r  =  g.  Then  by  the  definition  of 
B^*  (T),  we  have  g  €  B gm  (T)  with  3 *  satisfying  (2.2).  Then  by  Theorem  3.2  in  [4]  problem  (2.4)  has 
a  unique  solution  u  €  B^(D)  (up  to  a  constant)  with  3  satisfying  (2.6).  Then  applying  Theorem  2.1 
we  have  u|p  €  Bj  (T)  with  3  given  by  (2.1).  ^ 


Rc..tk  3.1.  In  general,  we  have  3  —  j  +  e  with  some  (  >  0.  According  to  (2.1).  i  2.2).  and  i  2.6 ,  J  (  . 

iresp.  Jl  2)  depends  on  the  interior  angle  and  on  j,  as  well  as  J,_,  iresp.  3^.).  For  instance,  if 
o  o 

g  6  B^i.r,)  n  0  <  3t,v  J,  have  —  max(  Jt  .j^j  -r-  ^  by  Stnre 

0  <  <  2tr.  i  -  4-  <  \  <  3'  and  J,  =  J*  by  (2.6V  Due  to  (2.1)  0,  ,  G  (J.-lJi.  hence 

l  ^  —  2  =  max( -  t-'i  k.  For  the  other  cases  the  relations  between  3t.  3 1.  J  _.  iresp. 

J1  +  1),  and  j-i,  can  be  derived  similarly. 

Next  we  derive  a  boundary  integral  equation  to  solve  the  Neumann  problem  1 2.1).  Insertinc 
the  fundamental  solution  v  =;  —  i  fn|x— y|  of  Au  =  0  into  the  second  Green  formula 

J (Auv  -  u Av)dx  =  J (u  p-  -  v  f^lds  ,3.1) 

Q  [" 

we  get  the  representation  formula  for  x  €  Q 


u(x) 


A  l  fn|x -y[u(y)ds(y ) 


J_  f  du(>~) 
2t  J.  dny 


£n|x  —  yjds(y ) 


which  yields  on  T  the  integral  equation 


(3.2) 


Du  =  (i_K')^ 

On 


(3.3) 


with  the  integral  operators  for  x  €  T 


and 


°U(X)  =  "  *  St  l  £->  ^|x-y|u(y)ds(y). 


du/  \  _  1  d  r  c  i  i  j  ,  \ 

K  ^(x)-  [  £nlx_yl  -q—  ds(y). 


Now,  equation  (3.3)  and  the  boundary  conditions  in  (2.4)  lead  to  the  First  kind  integral  equation  (see 

[12]) 


Du  =  f  on  T  (3.4) 

where  f  =  (1— K^)g.  There  holds  the  following  result: 

0,1 

Theorem  3.3.  For  given  g  €  (T)  satisfying  (2.5)  with  0  as  in  Theorem  2.4,  the  integral  equation 


(  3 . -4  )  together  with  the  side  condition 


J  u  ds  =  0  f  3  i 

r 

i.: 

has  a  unique  solution  u  £  B3  (T)  with  J  determined  by  (2.1)  and  (2.6). 

Proof.  Firs,  from  Remark  3.2  and  Lemma  3.2.  below,  we  have  g  £  H~1/2'T). 

By  Theorem  1.5  (ii)  of  [12]  together  with  uniqueness  of  the  system  (3.4).  (3.5)  .  its  solution  u 

exists  uniquely  in  H  (T).  Inserting  this  u  and  the  boundary  condition  g  —  ~  into  (3.2  i  we  define 

on 

for  x  £  Q 

Qfx)  =  'Jr  f.  chil  Cn\x-rMy)ds(y)  -  ^  J  g(y)£n|x- y!ds(y).  (3. 61 

By  il2]  u  £  H ‘ ( Q )  . solves  the  Neumann  problem  (2.1).  Due  to  the  derivation  above  v  =  u'p  £ 

1/2. 

H  ~ (  r )  solves  the  integral  equation  (3.4). 

0.1 

On  the  other  hand  from  Theorem  3.2  we  know  that  for  g  6  B-,  (T)  there  exists  a  unique 

solution  (up  to  a  constant)  l  £  B j ( Q )  of  the  Neumann  problem  (2.4)  with  J  satisfying  (2.6:.  But 

from  the  definition  of  B j( Q )  we  have  U  £  H*(f;)-  Hence  0  -  U  =  constant  on  Q.  But  for  v  =  u  r 

we  have  u  —  v  =  constant  on  T  where  u  £  H‘/‘(r)  solves  the  integral  equation  (3.4).  Hence 

...  1.2 

u  —  I  p  +  c.  c  —  constant,  and  since  l|p  £  Bg  (T)  due  to  Theorem  3.2  for  J  satisfying  (2.1),  the 
1.2  1 

assertion  u  £  B3  (T)  follows.  □ 

In  the  proof  of  Theorem  3.3  we  made  use  of  the  following  results. 

Lemma  3.1.  Lei  g  £  H=  (I)  with  0  <  3  <  5,  then  g  £  (H1/2(I))'. 

Proof.  We  may  assume  that  I  =  [0,1],  and  <j>;  =  x^1,  x  e  I,  0  <  /?,  <  5.  For  any  v  £  H 1 7 2 ( I ) 


Since  0  <  3l  <  |,  there  exist  integers  p,q  such  that  i  +  i  =  1,  and  2tf,p  <  1.  Therefore  by  the 
imbedding  theorem  (see,  e.g.,  [16]) 

K’Nvitix  <  (|  l-,<»’dX)‘/-(||V|!'dx),/'  <  C||V|£J<(I)  <  ci:v|iH,/I(i) 


9 


and 


ijgvdxj  <  Cilgfi  o  o  !|vji  .  /2 

I  I)  H  (I) 


which  leads  to  the  conclusion. 


Lemma  3.2.  7/g  E  H ]  iT).  0  <  d  <  then  g  G  H~i/2(F/. 


Proof.  First  from  [26]  we  observe  that  for  v  g  H1/2(r)  there  holds 


M 

£  IMi  1/2  <  ci|V|| 

i=!  H1/2(r,)  -  h  /_( r ) 


j 

with  a  constant  C  independent  of  v  and  the  length  of  Tt.  Hence  we  obtain  for  g  G  H  t  t  T  > 

Jgvds 


6l  H"1/2(r) '  IMI  1/2 


J  gvds 

r  Mr  M 

— -  <  C  £ 


<  C  £  lig||  o.o  <  C  g 


n  0.0 


h1/2(D  i  =  1"  l!H1/2(r,)  i  =  1  Rb  <r.)  h-  tD 


with  a  constant  C.  Here  we  have  made  use  of  Lemma  3.1  and  taken  f,  =  I.  Hence  by  the  definition 


of  H  l/‘(r)  the  assertion  follows. 


□ 


o,i  o 

Remark  3.2.  Definition  shows  that  g  E  B^  (T)  belongs  to  BDF)  with 


3, 


9* 

'3t~\ 


for  0  <  d,  <  | 


for  \  <  j,  <  1 


Furthermore  we  have  B^(T)  C  H£  (T)  due  to  the  definition  of  B|(T)  and  therefore  by  Lemma  3.2  we 
have  g  6  Bj  (T)  belongs  to  H_1/J(r). 


Now  we  discuss  the  numerical  solution  of  the  system  (3.4),  (3.5)  by  the  h-p  version  of  the 
Galerkin  boundary  element  method. 


Let  fi  be  an  L-shapea  domain  as  shown  in  Fig.  3.1.  We  assume  for  simplicity  that  the  solution 

>b  =  lyl  '  , 

02 


u  of  (3.4),  (3.5)  belongs  to  B^T)  (resp.  B^(T))  with  4>-  =  |y|^12,  <j>-  -  |x|5,'\  |  <  0U2.07A  < 


(resp.  |  <  0,  2,  /?2ii  <  1)  and  =  1  for  3  <  j  <  6, 
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1  1 


Fig.  3.1.  L-Shaped  Domain  Q. 

1/3  P 

i.e..  the  singularity  occurs  only  at  the  origin.  For  example  this  is  the  case  of  u  =  r  sin  ^  on  T  iresp. 

o  /3  O 

u  =  r*  cos  j  9)  where  (r,9)  denote  the  polar  coordinates  centered  at  the  origin. 


Let  o  €  (0,1)  be  the  mesh  factor  and  n,  integer,  be  the  number  of  layers,  and  let  T  ,  1  <  i  < 
I(j),  1  <  j  <  n  +  1  be  the  boundary  elements  such  that  dist(0,r,  )  -  an  +  ]  ;.  1  <  j  <  n  +  1  and 
dist(0,r,  j )  -  0.  1  <  i  <  l(j).  Then  f"  =  {T, 1  <  i  <  I(j),  1  <  j  <  n  +  1}  is  called  the  geometric 


11  11 
Mesh  ft  Mesh  tI 
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n  =  3 


1  1 
Mesh  r* 


Fig.  3.2  Geometric  Mesh  Fff,  n  =  1.2.3,  a  —  0.15. 


n=3 


Fig.  3.3.  Geometric  Mesh  n  =  3.  a  —  0.15 

Let  P  =  {p, 1  <  i  <  I(j),  1  <  j  <  n-4-1}  be  the  degree  vector  with  p,  ;  >  1  integer.  The 
boundary  element  space  associated  with  the  geometric  mesh  Y'l  and  degree  vector  P  is  defined  as 

Sp(r")  =  {d>  j^|p  is  a  polynomial  of  degree  <  p,  •  } 

and 

sp(r")  =  sp(r")  n  c°(r)  c  h,/2(d. 

The  geometric  mesh  r„  can  be  extended  to  a  geometric  mesh  Q"  on  as  shown  in  Fig.  3-3. 
Therefore  the  geoemtric  mesh  f"  on  dfi  can  be  defined  as  the  trace  of  the  geometric  mesh 


i.e..  f,  ;  =  50,  }  fl  <90.  Hence  the  space  SP(  )  is  actually  the  trace  spaces  of  Spf 0,  j.  i.e.. 


Sp(r")  =  {o!r:  0  £  Sp(0")} 

and 

sp(r")  =  {o|r.-  c  e  sp(o")  n  c°(0)}. 

For  details  of  the  definition  of  the  geometric  mesh  on  a  curvilinear  polygonal  domain  0  and  the  finite 
element  space  SP(0„),  see  [5],  [7],  [17]. 

From  [17]  we  quote  the  following  approximation  property  of  SP(Q"'. 

O  a 

Lemma  3.1.  Let  u  £  B^(0)  with  <J >5  =  r  .  0  <  3  <  1.  For  any  a  £  (0.1)  there  exists  w  £  Sp(0„) 
H  C°(0)  with  p,  j  —'Pj  >  1  and  jp  <  p;  <  nu.  0  <  p  <  v  <  oo,  which  satisfies 

""-|!B'(n, s  Cc_‘" 

where  n  is  the  number  of  layers  of  the  geometric  mesh  0„,  and  C,b  are  some  constants  depending  on  J 
and  a  but  not  on  n. 

By  the  trace  theorem  we  know  that  SP(T„)  C  H1/2(T).  The  corresponding  boundary  element 

■  01 
Galerkin  procedure  for  the  integral  equation  (3.4),  (3.5)  reads:  For  given  g  £  B;,  (T)  find 

Up  €  Sp(r")  such  that 

(Du„w)l2  -  ((l-K'jg.w^  V  w  £  SP(r^)  (3.7) 

and 

J  Upds  =  0 

r 

—1/2  1/2 

where  (•,•)  denotes  the  duality  between  H  (T)  and  H  (T).  We  have  the  following  approximation 
result  of  the  boundary  element  Galerkin  method. 

1  2 

Theorem  3.4.  Let  u  €  B^(T)  (resp.  B^(T))  be  the  solution  of  the  integral  equation  (3.4),  and  T  be  the 
boundary  of  the  L-shaped  domain  as  shown  m  Fig.  3.1,  where  0 ij  —  0xl  =  /?2  2  =  0,  3  <  t  <  6, 
J  =  1,2  and  g  <  012,  0 2,1  <  }  (resp.  |  <  0l2,  02, 1  <  !)•  Let  T*,  a  £  (0,1)  be  the  geometric  mesh 
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on  P:  and  let  SP(r,j)  denote  the  boundary  element  space  defined  above  with  p,  =  p,  >  1, 
i"  <  p;  <  nu,  0  <  /i  <  v  <  oc.  Then  the  boundary  element  Galerkm  solution  up  of  (3.7)  converges 
to  u  in  H  ‘(F)  exponentially,  i.e.. 


|u  —  Up| |  ,,,  <  Ce 

H  '  (f) 


-bN 


.1/2 


(3.8) 


where  N  is  the  number  of  degrees  of  freedom,  C  and  b  are  some  constants  depending  on  a  and  3  but 
not  on  N. 


1/2  • 

Proof.  Note  the  operator  D  in  (3.4)  is  strongly  elliptic  in  H  ‘(T),  i.e.,  satisfies  a  Garding  inequality 

in  H  ‘(T)  (cf.  Theorem  1.15  (ii)  in  [12]).  Hence  due  to  [19],  [27]  any  conforming  Galerkin  scheme 
1/2 

converges  in  H  (T),  and  we  have 


u  — u 


P  HI/2(f) 


<  C 


in;  |  u  Wn  1/2 

p,rV  PIIH  /2(f) 


wpesF(r 


Let  U  €  B^(G)  be  the  solution  of  the  boundary  value  problem  (2.4).  (2.5),  be  the 
geometric  mesh  on  Q.  ard  Sp(fit7)  fl  C°(fi)  be  the  2D  finite  element  space.  By  Lemma  3.1  there  exists 
a  Wp  €  Sp(tt")  n  C°(ft)  such  that 


|,u-wpMh1(Q)  <  Cie_6in  <3-9) 

where  n  is  the  number  of  layers  of 

Let  wp  =  Wpjr,  then  wp  e  Sp(r"),  and  by  the  trace  theorem 

l|u-wP|lHl/2(r)  =  |[U|r  -  Wp|rllHl/2(r)  <  C||U-Wp]|Hl(Q).  (3.10) 

Now,  (3.9)  and  (3.10)  together  yield 


-w. 


H1/2(r) 


<  Ce 


-6i  n 


Note  that  I(j)  <  k  uniformly  with  respect  to  j,  1  <  j  <  n+1.  Hence  N  <  kn2,  where  k  is  some 
constant  independent  of  N,  which  leads  to  (3.8)  immediately.  0 
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4.  Boundary  Element  Method  for  the  Dirichlet  Problem 

In  this  section  we  consider  the  Dirichlet  problem  (2.3). 

Lemma  4.1.  tf  g  E  B^  (T),  then  g  E  H1/*(I'). 

Proof.  Since  g  E  Bj  (T).  g|r  €  H3  ( T, )  for  0  <  j,  <  \  o-  g'r-  €  H 3  ( T, )  for  4  <  j,  <  1. 

1  <  i  <  M.  By  Theorem  2.2  there  exists  a  G  €  with  3‘  satisfying  (2.2)  such  that  G  j  =  e 

By  the  definition  of  Hj*  (Q),  G  €  H1(Q).  Therefore  g  E  H 1  / ‘ ( T ).  □ 

1.2  2 
Theorem  4.1.  Let  g  €  B3  (P)  Pi  C°(r).  Then  there  exists  a  unique  solution  u  E  B^( Q)  of  { 2.3)  with 

3  satisfying  (2.6).  Further  G  B3  (T)  with  3  satisfying  (2.1). 

anlp 

2 

Proof.  By  Theorem  2.3  the  problem  (2.3)  has  a  unique  solution  u  G  B^Q)  with  3  satisfying  (2.6). 

Therefore,  D°u  G  B^G)  for  [a|  =  1.  Applying  Theorem  2.1  we  obtain  G  B3  ( T )  with  3 

c/nip 

satisfying  (2.1).  D 


Next  we  derive  a  boundary  integral  equation  for  (2.3).  We  introduce  the  integral  operator 

V  §J(x)  =  -  I  J  ^7  <n|x-y|ds(y),  x  G  T  (4.1) 

and 

Ku(x)  =  -  |  J  £n|x  — yju(y)ds(y),  x  G  T  (4.2) 

P  0 Tly 

Taking  the  limit  x  G  T  in  (3.2)  we  obtain  via  the  well-known  jump  relation  for  the  double  layer 
potential  K  the  integral  equation  on  f 


V  |H(x)  =  (1+K)u(x).  x  €  T. 


Insertion  of  the  boundary  condition  of  (2.3)  into  (4.3)  leads  to  the  first  kind  integral  equation 


V  =  f  on  T 

an 


with  f  =  (1  +K)g  for  which  there  holds  the  following  result. 
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Theorem  4.2.  Lei  cap(F)  ^  1  where  cap(T)  is  the  capacity  [or  conformal  radius)  of  T.  Then  for 
given  g  E  Bi  (T)  O  C°(r)  there  ezists  exactly  one  solution  £  B3  (T)  of  the  integral  equation  M.li 

with  J  given  by  (2.1)  and  (2.6). 


Proof.  From  Lemma  4.1  we  have  that  g  6  Bj  (D  implies  g  6  H  (D.  By  Theorems  3.5.  3.8.  and 


3.9  in  [13]  the  integral  equation  (4.4)  has  a  unique  solution  £  H_/‘ 

1/"  ri 

g  6  H‘  *(F).  Then  inserting  this  together  with  the  given  data  g  =  u  into  (3.2)  we  define 


( F )  for  given  data 


a(x)  =  ^  J  ^n|x— y ]g(y)ds(y)  -  ^  J  |'(y)£n|x-y|ds(y),  V  x  €  Q. 


(4.5 ) 


By  Theorem  3.9  in  [13]  we  know  u  £  H!(fi)  and  solves  the  Dirichlet  boundary  value  problem 


(2.3).  Also  p 

om 


—  1/2 

€  H  '(F)  satisfies  the  integral  equation  (4.4)  due  to  the  definition  of  u.  Due  to  the 


uniqueness  of  the  solution  of  (4.4)  in  H 


-1/ 


(F)  havefS  =  i- 


and 


d\J\ 

<9n 


By  Theorem  4.1  there  exists  a  unique  solution  U(x)  6  By(Q)  of  (2.3)  with  3  satisfying  (2.6). 


0,1 


€  (F)  for  0  satisfying  (2.1).  Since  Bp(fl)  C  H’fQ)  we  have  u  =  U  by  the  uniqueness  of 


the  solution  of  the  Dirichlet  problem.  Hence 

on' 

completing  the  theorem. 


-  m 

dn\ 


0.1 


£  B5  (T).  But  we  also  have 


_  du 
dm 


Now  we  consider  the  rate  of  convergence  for  the  h-p  version  of  the  boundary  element  Galerkin 
method  for  the  integral  equation  (4.4).  For  simplicity  we  assume  again  that  Q  is  the  L-shaped  domain 
shown  in  Fig.  3.1  with  cap(T)  ^  1  and  that  g  and  the  solution  ^  of  (4.4)  have  a  singularity  at  the 
origin  only.  Then  the  geometric  mesh  on  T  and  the  boundary  element  space  SP(r„)  are  defined  as 
in  the  previous  section.  They  are  the  traces  of  the  geometric  mesh  Q„  on  Q  and  of  the  finite  element 
space  SP(Q„),  respectively.  Obviously  Sp-l(r")  C  L2(T)  C  H _  1  / 2( T). 

1.2 

The  Galerkin  procedure  for  the  integral  equation  (4.4)  reads:  For  given  g  £  Bj  (T)  n  C°(r) 
find  t lip  £  Sp_1(r")  such  that  for  all  4>v  £  SP-1(Fff) 

(Vtpp'<l,p)L\r)  =  ((1  +  K)s^p)L2(r)-  (46) 

For  the  boundary  element  solution  ipp  we  have  the  following  approximation  theorem. 

Theorem  4.3.  Let  gjj  €  Bj(T)  ( resp .  B^(T))  he  the  solution  of  the  integral  equation  (4.6)  where  T  with 
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cap(T)  ^  1  t s  the  bon  dary  of  the  L- shaped  domain  as  shown  in  Fig.  3.1.  and  2  =  J,  ,. 

3  <  i  <  6.  j  =  1,2.  1  <  j12.  3,,  <  ;  ( resp .  |  <  jl2.  J2l  <  1).  Let  a  €  (0,1)  be  the 
geometric  mesh  on  T  and  SP  1  ( r<, )  be  the  boundary  element  space  defined  in  previous  sections  with 
pt  =  p;  >  1.  jpi  <  p;  <  i/n,  0  <  p  <  i/  <  oc.  Then  the  boundary  element  Galerkin  solution  v P  of 
(4.6)  converges  to  ^  in  H  ^‘(P)  exponentially ,  i.e.. 


p 


9u|  <  Ce-bN'1/2 


(4.7) 


where  N  is  the  number  of  degrees  of  freedom,  and  C,  b  are  some  constants  depending  on  ty  and  J,  but 
not  on  N. 

Proof.  Since  the  operator  V  is  strongly  elliptic,  i.e.,  satisfies  a  Garding  inequality  in  H  (T)  (cf.. 
Theorem  2.19  in  [13]),  any  conforming  Galerkin  scheme  converges,  and  we  have 


3u  _ 

<9n  Vp 


H'1/2(r) 


<  C 


inf 


wpesp_1(r")1 


du 

dn 


-  w. 


H 


-1/2 


(D 


(4-8) 


Let  V  be  the  harmonic  conjugate  of  the  solution  U  of  the  boundary  value  problem  (2.3).  Then 
V  satisfies 


AV  =  0  in  fi 


dV\ 

ds\ 


r 


(4-9) 


and  V  €  B^(fl).  Following  the  proof  of  Theorem  5.1  of  [7]  we  can  show  that  the  2D  finite  element 
solution  Vp  €  Sp(fi")  H  C°(D)  of  the  Galerkin  equations  for  (4.9)  satisfies 


4,  n 


(4.10) 


where  C1,b1  are  independent  of  n  and  p.  Furthermore 


Vp  ds 


9V'I  =  k,  6  s'-‘(r:), 
r 


where  kp  is  the  projection  of  k  on  Sp  1(r<r).  Therefore  we  have  with  ^ 


(compare  Theorem 


4 . 2 '  and  using  (4.9) 


Idu  v  1 

<9U|  _  d\  P 

(dX 

dX  P\ 

bn  PIh_1/2(D  " 

<9njp  ds 

r 

H"i/2(D 

Us 

ds  ) 

r 

H~1/2 

From  Lemma  3.2  in  [25]  it  follows  that  w  £  H 1  /_ ( T )  implies  for  the  derivative  w'  £  H-l/'(n.  Hence 
there  exists  a  constant  C  such  that 


$¥  -  IT2 1  -1/2  <  C||V  — Vp||  1/2 

ds  ds  ly  1/2{D  _  H v  (T) 


Therefore  together  with  the  trace  theorem  we  obtain  from  (4.10)  and  (4.11) 


—  -  v  J 

dn  1/2(D  ~ 


<  C||V-Vp|!  .  <  Cle'‘in. 


(4.12 


Now.  (4.12)  together  with  (4.8)  yields  (4.7)  by  noting  that  N  <  kn2  for  some  constant  k  independent  of 

N.  □ 


Remark  1.  The  regularity  of  the  solutions  of  the  boundary  integral  equations  for  mixed  boundary 
value  problems  and  the  exponential  rate  of  convergence  for  the  n-p  version  of  the  boundary  element 
Galerkin  method  can  be  proven  similarly,  but  we  will  not  elaborate  it  here,  see  [4],  [6],  [5],  [7],  [13]. 

Remark  2.  In  previous  sections  we  assumed  that  Q  is  a  straight-line  polygon.  If  Q  is  a  curvilinear 

polygon  with  a  piecewise  analytic  boundary,  the  solutions  of  the  boundary  value  problems  (2.3)  and 
2 

(2.4)  belong  to  B^+<(Q)  for  any  e  >  0,  (see  [5]).  Then  it  could  be  proven  that  the  solutions  of 

1,2  o.i 

equations  (3.4)  and  (4.4)  are  in  B^+{(T)  and  B^+i(r)  respectively.  Thus  Theorems  3.4  and  4.3  remain 
valid. 


Remark  3.  The  geometric  mesh  shown  in  Fig.  3.2  is  designed  for  the  problems  with  a  singularity  at 
only  one  corner.  But  it  is  not  difficult  to  generalize  the  above  technique  to  the  problem  with 
singularities  at  each  corner  of  T,  and  the  exponential  rate  of  convergence  of  boundary  element  Galerkin 
solution  can  be  proven  again. 
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